Atoms and molecules in lattices: condensates built on a shared vacuum 
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In optical lattices where each site is occupied in its lowest energy state by a superposition of 
zero, one and two atoms, one can in a controllable manner convert the atomic pair into a molecule 
while retaining the vacuum and one-atom amplitudes. The microscopic quantum coherence on each 
site between the vacuum and the single molecule component leads to a macroscopically populated 
molecular condensate when the lattice is removed. 
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In this paper we present a simple quantum optical anal- 
ysis of the diffraction pattern of atoms and molecules re- 
leased from a periodic potential. We show that removal 
of atom pairs from a lattice with low filling fraction may 
be used to reestablish a robust atomic interference pat- 
tern after dephasing has taken place due to collisional 
interaction. If the removal occurs by photoassociation, 
a molecular condensate may be produced and detected 
when the lattice is removed. We do not make use of the 
Mott-insulator phase transition dynamics [Q, ^| . On the 
contrary, the existence of a superposition of vacuum and 
one or several atoms or molecules on every lattice site is 
crucial, and our proposal works ideally in the absense of 
atomic interactions. 

Assume that an atomic Bose-Einstein condensate has 
been exposed to a lattice potential, which was turned 
on slowly enough for the spatial motion of the atoms 
to remain in the lowest energy state but too fast for the 
Mott insulator transition to have happened. This system 
is described as a single condensate in the q = state 
of the lowest Bloch band of the potential, which is, in 
case of a deep potential modulation, a superposition of 
Wannier states localized in each well. We assume that 
all the atoms in the m'th well are well described by the 
m'th Wannier function <j) m {x) = <f>a(x — mL). 

If the entire many body wavefunction is described as 
a coherent state, i.e., an eigenstate of the atomic an- 
nihilation operators ^(x), the state remains a coher- 
ent state, which factors exactly into a product state 
of coherent states populating the Wannier states, i.e., 
each site is populated by a superposition *^2 n c n \n) m 
of zero, one, or more atoms, where the amplitudes c„ 
can be parametrized by a single coefficient f3: c n — 
cxp(— \[3\ 2 /2)f3 n /\fri\. Such superposition states do not 
reflect the conservation of atom number. The fluctua- 
tions in atom number on each lattice site, however, re- 
flects the entanglement of the wells due to the splitting 
of the atoms - and it is a remarkable feature of coher- 
ent states that this entanglement is well accounted for 
by a product state (because a formal projection on a 



total number eigenstate with very many atoms, which 
indeed produces a non-separable state, does not change 
our physical predictions). 

Since atoms are bound to occupy the ground states of 
the wells of the lattice, it is convenient to introduce the 
discrete set of atomic field operators that removes atoms 
from the Wannier mode functions rather than from a 
specific location in space. These arc defined as 
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and they obey the standard commutator relations 

With the product wave function = 
Cn|?^)m); the mean value of the mode anni- 
hilation operator in each well equals a := (a m ) = 
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number of atoms in each well 
= Y^=i\ c n\ 2 n, and we have 
the inter-well coherence (a^am) = a*a for m ^ m! . 
The short range interaction between atoms is a weak 
perturbation which causes a phase evolution of the c n 
amplitudes with a frequency ^-n(n — 1). In the absense 
of tunneling between wells, the product form of the wave 
function remains exact, but the phase evolution leads 
to a reduction of the atomic field amplitude a since the 
terms c^^Cn^fn acquire different complex phases. The 
resulting disappearance of atomic interference has been 
observed in experiments [^[ |J, and subsequent revivals 
when the phase differences reach multiples of 27r have 
also been shown pjj. 

We now turn to the special case of on the order of or 
less than unit mean occupancy of each well. The am- 
plitudes on states with three or more atoms are small 
and we shall mainly focus on the lower n components. 
Since \cqCi +c\c 2 V2| may be smaller than |cqCi|, it seems 
an interesting experiment to remove the component with 
two atoms from the state vector of each well. This can 
be done by a carefully designed photoassociation process 
P, 0, Q], making use of the fact that the two-atom com- 
ponent is in a fully determined initial state. Neglecting 
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higher atomic occupancies, the resulting state vector be- 
comes a product of two-species superposition states: 

|*) = n m (c |vacuum) + Ci|l atom) + c 2 |l molecule)) m . 

(2) 

In the state (^|) we can readily determine the mean 
atomic amplitude and population, a — CqCi and n — 
\c\\ 2 . The mean amplitude acquires a finite robust value, 
which is not corrupted by the atom-atom interaction be- 
fore or after the photoassociation. We may thus wait 
arbitrarily long before we remove the atomic pairs, and 
still the coherence reappears. Technically, we define con- 
densate population as the largest eigenvalue of the one- 
body density matrix, and applying the basis of Wannier 
states, this matrix has the mean occupancy of the wells 
n in the diagonal, and the squared amplitude \a\ 2 in all 
other positions. In a system with N wells, the eigen- 
vector with equal amplitude on each Wannier function 
has the largest eigenvalue of (N — 1)M 2 + n, i.e., in the 
limit of large N the condensate fraction is \a\ 2 jn. This 
is also the value predicted by the off-diagonal long-range 
order, (a^ami) = \a\ 2 = n ■ (\a\ 2 /n). It is interesting to 
note that in cases where c 2 is sufficiently dephased, pho- 
toassociation increases the value of \a\ 2 and hence both 
the total occupancy of the condensate and the condensate 
fraction increase by the removal of atoms. The state be- 
fore photoassociation is a perfect, pure quantum state of 
the system, and it is related to the perfect condensate 
state by a unitary process - but the system cannot be 
described by a Gross-Pitaevskii wave function, and the 
one-body density matrix cannot distinguish a determin- 
istically evolved phase from, e.g., a finite temperature 
decoherence effect. 

Let us now turn our attention to the molecules pre- 
pared in the lattice. It is clear from the state (|2|), that 
the atomic and molecular states are quite equivalent: 
there is a molecular population um = | C2 1 2 per lattice 
site, and there is a mean molecular field of au — CqC2- 
The above analysis of condensate population and frac- 
tion thus applies to the molecules if n and a are replaced 
by um and ccm- The preparation of molecular conden- 
sates has received quite some interest, in particular with 
the recent experiments at JILA || , where coherent oscil- 
lations between an atomic condensate and a condensate 
of very loosely bound molecules are observed and agree 
well with a detailed theory of collision dynamics near a 
Feschback resonance JlCfl , see also (yj. In addition to 
the detailed binary collision dynamics, it is necessary to 
represent correctly the many-body dynamics in the case 
where many atoms are converted into molecules. In com- 
parison, our analysis is much simpler since there is either 
nothing or just a single atom or just a single molecule in 
every lattice site (as long as the truncation of the atomic 
state above two atoms is a valid approximation). There 
is accordingly just a single population independent cou- 
pling between the relevant states \n = 2 atoms) and 



\um = 1 molecule). From an experimental point of 
view our molecules can be prepared in a deeper bound 
molecular state, selected coherently by the coupling laser 
fields rather than by Feshbach collision dynamics, and 
they are prepared in a spatial condensate for which the 
phase coherence can be readily detected by diffraction 
experiments similar to the ones applied to atoms in [0, 
see below. 

Our proposal may also be compared with the sug- 
gestion by Jaksch et al. where the Mott-insulator 
dynamics is used to prepare a lattice with precisely 
two atoms per site, these atoms are transferred into 
a molecule, and the molecular Mott-insulator is sub- 
sequently melted to yield the molecular condensate. 
Our proposal emphasizes a simple but spectacular phe- 
nomenon: a zero-quantum amplitude is enough to yield 
a mean field, and hence the condensate establishes itself 
simply by the production of molecules out of the two- 
atom component of the atomic state. This is so, because 
the atomic vacuum is of course also a molecular vacuum 
(as well as a vacuum for larger molecules, butterflies and 
freight trains). On the experimental side, the molecular 
condensate may ideally be created from an atomic con- 
densate of non-interacting particles, and we only need to 
transfer the atoms slowly enough into the lattice to have 
a well-specified initial state for the photo-association pro- 
cess. 

The analysis is quite simple and fully analytical in the 
limit of very weakly populated sites. Also for a larger 
mean number of atoms on each site the conversion into 
molecules will lead to a molecular mean field, i.e., a 
molecular condensate. To investigate this, we have solved 
numerically the time evolution of the simple photoasso- 
ciation Hamiltonian 

Hp A = X(a 2 m bl + (al) 2 b m ) (3) 

m 

where b m is the operator of annihilation of a molecule at 
site 77i. Starting with a coherent atomic state, the Hamil- 
tonian (||) introduces in every lattice well a superposition 
state Y^, c n,n M Wl n > n M), from which the mean number of 
molecules and the mean molecular field is readily calcu- 
lated. The results of such calculation are shown in Fig.l. 
It shows the number of molecules and the absolute square 
of the molecular field as functions of time: the conversion 
is oscillatory in time, and to optimize the coherent molec- 
ular component, one should stop the photoassociation 
process after a time interval depending on the number of 
atoms per lattice site. For simplicity we did not incorpo- 
rate atom-atom interactions or interactions between the 
atoms and molecules in these calculations since our main 
focus is on lattices with low population. 

In the above analysis we have shown that the on-site 
mean field amplitude and population mathematically de- 
termine the condensate fraction. These quantities are di- 
rectly measurable by the diffraction occuring when the 
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particles are released from the lattice. The detection of 
the diffraction pattern has precisely been the hall mark of 
coherence ||, [|, and coherence decay and revivals ||| in 
experiments, and we present here a brief theoretical de- 
scription, linking the diffraction pattern directly to our 
quantities n and a. 

If one introduces a complete and orthonormal set of 
wavefunctions {</> m ./}, where 4> m .o — <f) m are the Wan- 
nier functions applied above, and {4> m ,i} with I > are 
higher excited states centered on the m th lattice site, the 
commutator relations [&(x), W{x')\ — 5{x — x'), lead to 
5(m, with the mode annihilation op- 
erators defined through 

a m ,i = f <t> m ,i(x)$(x)dx. (4) 

We shall need expressions for the position dependent an- 
nihilation operators in terms of the mode functions, 

= }^a m ,i^ m ,i( x ) 

= flmfe(l) + Si>(x), (5) 

where in the second line we include only explicitly the 
Wannier state <j>^, where to enumerates the well which 
is most close to the position coordinate x, and where 
the 'noise-term' S^(x) represents all other modes, which 
are either unoccupied (I > 0) or which have no support 
at the location x, (to ^ to). This noise term does not 
contribute to atom counting signals. 

In experiments on trapped condensates, the spatial 
properties of the atomic cloud are detected by releasing 
the cloud and letting it expand. In case of the weakly 
filled lattices, one may neglect interactions during the 
expansion, and the final spatial distribution hence maps 
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FIG. 1: Incoherent and coherent molecular components pre- 
pared by photoassociation of atoms. The number of molecules 
(solid line) and the squared norm of the mean molecular field 
(dashed line) as a function of time is shown, starting at time 
zero with a coherent atomic state with unit mean number. 



the momentum distribution of the atoms in the lattice. 
We determine the momentum distribution by writing the 
atomic field operator that annihilates an atom with mo- 
mentum hk, in terms of the spatial field operators: 

4f(k) = -£= [ dxe lkx V{x) 
V 2ir J 

= J2ame imL <l>(k)+^(k), (6) 

771 

where <p(k) is the Fourier transfom of the Wannier func- 
tion cf>o(x), and where <5 x I'(fc) is a vacuum noise field with 
zero mean and zero intensity. Although the quantity to 
originally appeared as a complicated index, depending 
on x, the integral over all positions leads to a sum over 
all positive and negative integer values for to, i.e., it has 
become a conventional summation index. The number 
of atoms with momentum hk is given by the expectation 
value of the operator ^/^ (k)^(k): 

= \Hk)\ 2 E e t(m - m ' )kL (ala m ). (7) 

m,m' 

Thus 

(*t(fc)$(fc)) 

= rnrni 2 - M 2 ) + m 2 E s ( k ' ' (8) 

where N is the (large) number of lattice sites. The sum 
over q gives rise to a comb at lattice momenta q^f- with a 
modulation proportional to the square of the mean field 
amplitude. This comb sits on top of a flat background 
caused by the incoherent population of the wells, and the 
whole distribution is comprised within the width of the 
single well momentum distribution |(/j(/c)| 2 . 

A diffraction pattern in the molecular distribution is 
the clear signature of a molecular condensate. By count- 
ing the released atoms or the molecules on a position 
sensitive detector the values of n and a and thus the 
condensate fraction can by determined. 

In summary, we have presented a method which by 
photoassociative removal of atoms can at the same time 
return an atomic system to a state with a macroscopic 
population of a single quantum state and prepare a 
molecular system with a similar macroscopic population. 
Weak fields owe their mean amplitude to the existence 
of both vacuum and single quantum excitations, and the 
molecular condensate exists due to the molecular vac- 
uum component being already populated in the atomic 
condensate. 

Eq. (||) relates the diffraction pattern to the mean field 
amplitudes existing prior to the detection. Two indepen- 
dently prepared condensates also show interference |L2| . 
It has been demonstrated theoretically, jl3| , that this in- 
terference builds up as a consequence of the first random 
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detection events, and that one can think of the state of 
the two condensates without initial mean fields as a sta- 
tistical mixture of states with different relative phases 
Jilj , which all yield interference patterns but with differ- 
ent offsets. This same mechanism is not at work when 
atoms or molecules emerge from a large number of inde- 
pendent sites, since the uniform distribution of the phase 
variable at every site has vanishing probability to pro- 
duce the regular progression of the phase over several 
sites needed for a high visibility diffraction pattern. 

Our approach is quantum optical in nature: by care- 
fully selected unitary operations we generate from a given 
initial state a final desired state of the system. This im- 
plies that the formal properties of the state vector of 
the system are made to the design, but the dynamical 
properties of the system may have only little in common 
with the properties usually ascribed to a Bose-Einstein 
condensate. Since the state is not the ground state or 

the thprrnndynarnir pqnilihrinm state nnHnr an applipH 
Hanliilt.rmian it. will not act in the samp way as snrh a 
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